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Continuing our study on the complete integrability of nonlinear ordinary differ- 
ential equations, in this paper we consider the integrability of a system of coupled 
first order nonlinear ordinary differential equations (ODEs) of both autonomous 
■ and non-autonomous types. For this purpose, we modify the original Prelle-Singer 

procedure so as to apply it to both autonomous and non-autonomous systems of 
coupled first order ODEs. We briefly explain the method of finding integrals of 
motion (time independent as well as time dependent integrals) for two and three 
coupled first order ODEs by extending the Prelle-Singer(PS) method. From this 
we try to answer some of the open questions in the original PS method. We also 
identify integrable cases for the two dimensional Lotka-Volterra system and three- 
dimensional Rossler system as well as other examples including non- autonomous 
systems in a straightforward way using this procedure. Finally, we develop a lin- 
earization procedure for coupled first order ODEs. 

00 

O ' 1- Introduction 

00 '. 

In our previous two works (Chandrasekar et al. 2005; 2006) we have studied in some 
detail the extended modified Prelle-Singer (PS) procedure (Prelle & Singer 1983; 
Duarte et al. 2002) so as to apply it to a class of second and third order nonlinear 
ordinary differential equations (ODEs) and solved several physically interesting 
nonlinear systems and identified a number of important linearization procedures. 
We now wish to extend the procedure to coupled ODEs. In the present paper we 
discuss the modification and applicability of the extended PS method to a system 
of first order ODEs of both autonomous and non-autonomous types. In subsequent 
papers, we will extend the procedure to coupled second and higher order ODEs. We 
are motivated by certain open questions still prevailing in the original PS method 
for a system of autonomous coupled first order ODEs. Before discussing them and 
describing how we answer them, we shall have an overview of the original PS method 
and its generalizations. 

We recall that in 1983 Prelle & Singer (1983) have proposed an algorithmic 
procedure to find the integrating factor for the system of two autonomous first order 
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ODEs of the form ^ = P(x, y), = Q{x, y), where P and Q are polynomials in x 
and y with coefficients in the field of complex numbers. Equivalently it can be recast 
in the form y' = ^| — p|^j ■ Once the integrating factor for the latter equation 
is determined then it leads to a time independent integral of motion for the above 
autonomous two coupled first order ODEs (for the single first order ODEs the first 
integral is nothing but the general solution). The PS method guarantees that if 
the given system of two coupled first order ODEs has a first integral in terms of 
elementary functions then this first integral can be found. This method has been 
generalized to incorporate the integrals with non-elementary functions (Singer 1990; 
Singer 1992; Duarte et al. 2002), and some first integrals of autonomous systems 
of ODEs of higher dimensions (of dimension 3) were also calculated. Man (1994) 
described a method for calculating first integrals of autonomous systems which 
are rational or quasi-rational, but said that 'The generalization of this procedure 
to higher dimensions to find elementary first integrals is still an open problem'. 
In addition, the question whether the PS procedure can be extended to a non- 
autonomous system of first order ODEs has not been addressed so far. Also to our 
knowledge the problem of finding time dependent integrals for a given system of 
coupled first order ODEs through this procedure has not been taken up so far. 
Further, the problem of finding both time dependent and independent integrals 
for a system of first order ODEs greater than two in number has not been dealt 
with systematically. In addition, the problem of how to linearize a given system of 
coupled first order ODEs systematically also remains to be tackled. In this paper, 
we address positively all these questions and come out with valuable answers to 
these problems and demonstrate the general results with suitable examples. 

Firstly, in order to extend the PS method to a non-autonomous system of first 
order ODEs with the rational form £ = f = g^g, where P^,x,y) 

and Qi(t,x,y), i = 1,2, are analytic functions of x and y with coefficients in the 
field of complex numbers, we develop a modified technique and derive determining 
equations for the integrating factors R and K. In the case of a coupled system of two 
first order ODEs we observe that the integrals of motion ii and I 2 of the coupled 
ODEs are either both time dependent or I\ may be time independent while I 2 is 
time dependent (or vice versa). Based on this observation we analyze the problem 
by splitting it into two different cases, namely (i) the solutions of the determining 
equations for R and K which lead to time independent integral(s) and (ii) the 
solutions which lead to time dependent integrals. We show that the determining 
equations in the time independent integral case coincide with the determining equa- 
tions derived by Prelle and Singer in their original paper (Prelle & Singer 1983). In 
the second case we obtain the time dependent integral, if it exists, and from this we 
tackle the second problem. Using this method we find integrable cases for the two 
dimensional Lotka-Volterra (LV) equations (given as an example) with both time 
dependent and independent integrals as well as certain non-autonomous systems. 

For extending the method to a coupled system of more than two first order 
ODEs, we first extend the above modified PS procedure to three coupled first 
order ODEs and propose a systematic procedure to obtain the integrating factors 
R, K and M. In this case we divide our analysis into five categories. We illustrate 
this theory with two physically important examples, namely Rossler system and 
3D LV system and identify some new integrable cases. The method is extendable 
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straightforwardly to a system of more than three coupled first order ODEs. In 
addition to the above, we also present a method of finding linearizing transformation 
for a system of first order ODEs. We illustrate the theory with certain concrete 
examples. 

The paper is organized as follows. In the following section we develop the mod- 
ified Prcllc-Singcr method applicable for non-autonomous system of two coupled 
first order ODEs. In §3, we describe the method of solving the determining equa- 
tions and how one can obtain compatible solutions from them. In §4 we illustrate 
the procedure with LV equation as an example and identify many intcgrable cases 
in it. We also apply the method to non-autonomous two coupled equations. In §5, 
we extend the PS procedure to non-autonomous three coupled first order ODEs 
and describe methods of solving the determining equations in §6. We emphasize 
the validity of the theory with two illustrative examples arising in different areas of 
physics in §7. Further, in §8 we discuss the direct applicability of the modified PS 
procedure to n (>3) coupled first order ODEs. In §9, we demonstrate the method 
of identifying linearizing transformations with examples. Finally, we present our 
conclusions in §10. 

2. Prelle-Singer procedure for two coupled first order ODEs 

Let us consider a system of two coupled first order ODEs of the form 

Pi(t,x,y) P2(t,x,y) r - = j^ con 

Qx(t,x, y y y Q2(t,x, y y 1 dt ! [ ■ > 

where Pi and Qi, i = 1,2, are analytic functions of x and y with coefficients in 
the field of complex numbers. Further, we assume that the ODE (2.1) admits a 
first integral I(t, x, y) = C, with C constant on the solutions, so that the total 
differential becomes 

dl = I t dt + I x dx + Iydy = 0, (2.2) 

where subscript denotes partial differentiation with respect to that variable. Let us 
rewrite (2.1) in the form 

P P 

-±dt -dx = 0, -?-dt -dy = 0. (2.3) 
Qi Q2 

Hence, on the solutions, the f-forms (2.2) and (2.3) must be proportional. Multi- 
plying the first equation in (2.3) by the integrating factor R(t, x, y) and the second 
equation in (2.3) by a second integrating factor K(t,x,y) (both of which are to be 
determined), we have on the solutions that 

dI={R(f> 1 +K<p 2 )dt-Rdx-Kdy = 0, (2.4) 

where = Pi/Qi, i = 1,2. Comparing equations (2.4) and (2.2) we have, on the 
solutions, the relations 

It = {R<j)i + Kfa), I x = -R, I y = -K. (2.5) 
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The compatibility conditions, I tx — I x t, hy = lyt, Ixy — lyx, between the 
equations (2.5) provide us the conditions 

Rt + <t>iRx + foRy = -{R<j>i x +K<j> 2x ), (2.6) 

K t + <f>iK x + faK y = -{R<t> ly +K<j, 2y ), (2.7) 

Ry = K x . (2.8) 

Integrating equations (2.5), we obtain the integral of motion, 



I = r x +r 2 



-/[-!( 



n + r 2 



dy, (2.9) 



where 



n = fiRfa+K^dt, r 2 = -f(R+±(n))dx. 

Solving the determining equations (2.6)-(2.8) consistently, we can obtain expres- 
sions for the functions R and K. Substituting them into (2.9) and evaluating the 
integrals we can deduce the associated integral of motion. If two independent sets 
of solutions R and K for (2.6)-(2.8) are found then they give rise two independent 
integrals for the given system of first order ODEs (2.1) which ensures the complete 
integrability of the system and obtaining explicit general solution of the system. 

3. Method of Solving Determining Equations 

One may note that the determining equations (2.6)-(2.8) are over-determined and 
the crux of the problem lies in finding the explicit solutions satisfying all the three 
determining equations, since once a particular solution is known then the integral 
of motion can be readily constructed. To solve the equations (2.6) and (2.7) we 
introduce a transformation 

R = SK, (3.1) 

where S is a function of t,x and y so that the determining equations (2.6)-(2.8) 
now become 

s t + 4>\S X + 4> 2 Sy = -4>2x + {4>2 y - <t>ix)s + 4>i y s 2 , (3.2) 

K t + <f>iK x + <f> 2 K v = -K(Scp lv + <p 2y ), (3.3) 
K x = SK y + KS y . (3.4) 

One may observe that the equation for S, namely equation (3.2), is decoupled from 
that of K (equations (3.3) and (3.4)) and so the set (3.2)-(3.4) may be easier to 
analyze than solving the original ones (2.6)-(2.8) directly. 

To begin with we observe that for the system of two coupled first order ODEs 
(2.1) there can be two independent integrals I\ and I 2 such that either both of them 
are time dependent or I\ may be time independent while I 2 is time dependent (or 
vice versa) . So we consider these two cases separately with corresponding solutions S 
and K for equations (3.2)-(3.4). The determining equations in the time independent 
integral case (see equation (3.8) below) coincides with the determining equation 
derived by Prelle and Singer in their original paper (Prelle & Singer 1983), whereas 
for the time dependent integrals case we develop an extended procedure to capture 
both the integrals, if they exist, and thereby making the PS procedure a more 
powerful tool in a self-contained way. 
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(a) Time independent integral 

In the case I t = we denote I\= I and note that the function S can be easily 
fixed with the help of the first equation in (2.5), that is, 



R _ q _ 



(3.5) 



Since I is independent of t, it follows from equation (2.5) that S (and also R and 
K) is also independent of t. Indeed one can check that S — — 4^ is a solution of 

(3.2). Now substituting S = — ^ into (3.3) we get the following equation for K, 
that is, 

02 

(j>iK x + 4> 2 K y = K(—4> ly - <j> 2 y). (3.6) 
01 

We now make a substitution 

K = 7T-V ( 3J ) 

where f(x,y) is an arbitrary non-zero function of x and y. Then (3.6) takes the 
simpler form 

<t>lfx + fafy = f{4>lx+4>2y)- (3.8) 

We note that by redefining / = equation (3.8) coincides with the determining 
equation derived by Prclle & Singer (1983) for the autonomous case of (2.1). We 
also mention that equation (3.8) is nothing but the one obtained by substituting 
the forms of S and K given in (3.5) and (3.7) into the constraint equation (3.4). So 
by solving equation (3.8) we can get the complete set R(= SK) and K associated 
with equations (3.2)-(3.4). 

Even though equation (3.8) is a quasilinear PDE in two variables, the associated 
characteristic equation again leads to coupled differential equations of the form 
(2.1). Thus the routine methods of finding general solution of quasilinear PDEs are 
not very useful here. Here, we find particular solutions for the determining equation 
(3.8) in a different way and obtain integrable cases for the given system. For this 
purpose, we assume specific functional form for f(x,y) with unknown functions 
and determine the later consistently. A simple but effective choice is / = (A(x) + 
B(x)y) r , where A and B are functions of their arguments, and r is a constant. Again 
the reason for choosing this form is as follows. Since K is in a rational form, while 
taking differentiation or integration the form of the denominator remains the same 
but the power of the denominator decreases or increases by a unit order from that 
of the initial one. So instead of considering / to be of the form / = A(x) + B(x)y, 
one may consider a more general form / = (A(x) + B(x)y) r , where r is a constant 
to be determined. Depending on the problem in hand, one can also assume more 
general form and proceed as in the present case. 

Substituting now the form / = (A(x) + B(x)y)' r into (3.8) we arrive at the 
following equation for the unknown functions A and B, that is, 



h{A x + B x y) + <f> 2 B 



(<j>i x + foy)(A + By). (3.9) 
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Inserting the given form of 0,'s, i — 1,2, into the above equation (3.9) and solving 
the resultant equation, one can fix the forms of A, B and r. Now plugging the 
resultant form of / into equation (3.7) one can get the integrating factor K which 
in turn leads us to the other integrating factor R through the relation (3.5). Finally, 
substituting R and K into equation (2.9) and evaluating the integrals one can 
deduce the time independent integral for the given system. Since we are dealing 
a system of two first order ODEs, this time independent integral itself guarantees 
the integrability of the given system. However, to explore the general solution one 
may seek the time dependent second integral. We describe the procedure in the 
following sub-section. 

(b) Time dependent integral 

Now we focus our attention on the case 1% ^ 0. In this case, the function S has 
to be determined from equation (3.2). Since it is too difficult to solve equation (3.2) 
for its general solution, we seek particular solutions for S, which is sufficient for our 
purpose. In particular, we seek a simple rational expression for S in the form 

s= Mt,x)+B 1 (t,x)y 

A 2 (t,x) + B 2 (t,x)y 1 v ; 

where A^s and Bj's, i = 1,2, are arbitrary functions of t and x which are to be 
determined. Of course this can be further generalized, if need arises. Substituting 
(3.10) into (3.2) and equating the coefficients of different powers of y to zero, we 
get a set of determining equations for the functions A^s and Si's, i = 1,2. Solving 
these determining equations we obtain explicit expressions for the functions Aj's 
and Bi's, i = 1,2, which in turn fixes S through the relation (3.10). 

Now substituting the forms of S into equation (3.3) and solving the resultant 
equation one can obtain the corresponding forms of K. To solve the determining 
equation for K we again seek the same form of ansatz (3.7) but with explicit t 
dependence on the coefficient functions, that is, K = ^ t x ^g^ t x ^ y y , where Sd 
is the denominator of S. Once S and K are determined then one has to verify 
the compatibility of this set (S,K) with the extra constraint equation (3.4). Now 
substituting i?'s (= SK) and K's into equation (2.9) and evaluating the integrals 
one can construct the associated integrals of motion. 

We note here that for a given equation (2.1) one may also get two time dependent 
integrals or one time dependent and one time independent independent integrals 
(discussed earlier) which in turn automatically guarantees the complete integrability 
of the given system and provide us explicit solution by algebraic manipulation. On 
the other hand, under certain circumstances, one may get only one time dependent 
integral and one can transform this time dependent integral into time independent 
one and thereby establish the integrability. 

4. Two coupled ODEs - Application 

(a) Example: Two dimensional Lotka-Volterra (LV) system 

Our motivation is to identify integrable cases and deduce both time depen- 
dent and independent integrals for a given system of two coupled first order ODEs 
through the extended Prelle-Singer procedure in a self-contained way. 



Article submitted to Royal Society 



Integrability and Linearization 



7 



To demonstrate this we consider the celebrated two dimensional LV system 

x = x(ai + bux + bi2y) = <f>i, y = y(a 2 + b 2 ix + b 22 y) = 4> 2 , (4.1) 

where a^s and 6y's i,j — 1,2, are six real parameters. This system was originally 
introduced by Lotka (Lotka 1920) and Volterra (Volterra 1931) to model two species 
competition. However, in recent years this model appears widely in applied mathe- 
matics and in a large variety of physics topics such as laser physics, plasma physics, 
convective instabilities, neural networks, etc. (Minorsky 1962; Brenig 1988; Murray 
1989). The integrability properties of the system (4.1) alone have been analyzed by 
many authors, see for example Cairo & Llibre (2000) and references therein. 
In the following, we identify intcgrable cases in (4.1), through our procedure. 



(i) Time independent integral (I t = 0) 

In this case the function S can be fixed easily in the form (vide equation (3.5)) 

s= <t>2 = y(a2 + b 21 x + b 22 y) , 42 x 
4>i x(ai+bnx + bi 2 y)' 

To explore the integrating factor K we need to fix the form / first, see equation 
(3.7). For this purpose we substitute the (pi's, i = 1,2, into (3.9) so that one gets 
the following equation for the unknown functions A and B which constitute the 
function /, that is, 

x(a 1 + b n x + b 12 y)(A x + B x y) + y(a 2 + b 21 x + b 22 y)B 

= (oi +a 2 + (621 + 2b u )x + (612 + 2b 22 )y)(A + By). (4.3) 

Equating the coefficients of various powers oiy l ,i = 0, 1, 2, and solving the resultant 
differential equations for A and B we arrive at the following two general expressions 
which involve the system parameters, that is, 

b 2 \{b\ 2 - b 22 )(bi 2 b 2 i - 611(612 - 2622)) (ai (611 - 621)622 + »26ii(6 22 - 612)) = 0, 
011622(611 - 6 2 i) 2 (a 2 6i2 - ai6 22 ) + ^611(612 - b 22 )(b 12 b 2 i - 611622) = 0, (4.4) 

and 

r _ (611612 + 6 i2 6 2 i - 26n6 22 ) ^ _ (&11&12 + 6 i2 6 2 i - 2611622) ^ ^ 

611(612 - 622) (612621 - 611622) 

Any consistent solution which comes out from the above expression (4.4) gives us 
an integrating factor which in turn leads us to an integral of motion. In this sense 
(4.4) forms an integrability condition of some generality which in fact encompasses 
all known integrable cases with time independent integrals (Cairo & Llibre 2000; 
Llibre & Vails 2007). For example, while analysing the above equation, we find that 
one can straightforwardly recover several known integrable cases like (i) 6n = 6 22 = 
0,r = 1, (ii) b 22 = 612,611 = 621, r = 1, (Hi) a x = a 2 ,b 12 = 36 22 ,6n = -621, r = 
— 1, (iv) a\ = a 2 = 0, r as given in equation (4.5), and so on straightforwardly from 
equations (4.4) and (4.5) and construct the associated integral of motion which in 
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turn coincide with the existing results. However, as we are interested to construct 
an integral of motion with more general parametric choice we do not fix any relation 
between parameters (other then the general relation) and proceed further. 
The respective forms of A and B are (for b 22 bi 2 (b 22 — 612) 7^ 0), 



A=((r- l)g(x) + rb 22 (a 1 + b n x))x 
so that 

((2-r)t>22+i-12) 



((2-r)i, 22 +bi 2 ) 



B 



((2-r)b 22 +b 12 ) 

rb 2 2b 12 x ^ (4.6) 



where 



g(x) = a 2 b 12 - 2a 1 b 22 + 



(r - l)g{x) + rb 22 (a 1 + b n x + b 12 y) 



■(612621 + 611(612 - 26 22 ))a; 



(622 - 612 



(4.7) 



(4- 



Making use of the explicit forms of S and /, vide equations (4.2) and (4.7) 
respectively, with the parametric restrictions (4.4), the integrating factors, K and 
R, can be fixed as R = — -y- and K = y-. Substituting the forms R and K into (2.9) 
and evaluating the integrals we arrive at the following time independent integrals 
for (4.1) for the parametric cases (4.4) for different values of r, namely 



622 



<x) s 



622(01 + bux + 6i22/)(a 2 6i2 - 01622 



(612 - 622) 
26 22 (ai + b n x + b 12 y) 



(611622 - 612621)2; + 6126222/)) 



r^0,2, 



b 12 h(x,y) 



- log 



~h{x,y) 



I = x »i2 y(2a 1 + 2b 11 x + b 12 y) 1 



0, 



(4.9) 

(4.10) 
(4.11) 



where h{x,y) = a 2 + (6n + 621)2; + 2b 22 y and g{x) is given in (4.8). 

For 622612(622 — 612) = 0, one obtains known integrable cases following the 
same procedure as above. For example, let us consider the first case, which we 
cited above as the known case 622 = 611 = 0. For this case, one can find a trivial 
solution for equation (4.3) as A = and B = x with r = 1 so that / becomes 
/ = xy. The respective integrating factors R and K read as R = — a2 +^ 2lX and 
K = ° 1+bl2y so that the associated integral of motion takes the form / = 6 2 ix — 
bi 2 y + a 2 \ogx — cii logy. This integral is well known and popular in the literature for 
a long time (see for example Minorsky 1962; Prelle & Singer 1983; Murray 1989). 
Similarly for the case 622 = 612 and 6n = 621, we obtain the integral of motion of 
the form / = x a2 y~ ai {a\a 2 + a\b 22 y + a 2 bnx)( ai ~ a2 \ In the case 612 = 0, leads to 
the uncoupled equation for LV and the general solution for this case can easily be 
obtained. 



(ii) Time dependent integrals (I t =/= 0) 

Now let us concentrate on the case I t ^ 0. In this case S has to be determined 
from equation (3.2), that is, 

S t + x(a 1 + b n x + b 12 y)S x + y(a 2 + b 21 x + b 22 y)S y = -b 21 y 

+ ((a 2 + 6212; + 26 22 y) - (ai + 2b n x + b 12 y))S + b 12 xS 2 . (4.12) 
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As we mentioned earlier, to obtain a particular solution for the above equation (4.12) 
we seek a simple ansatz for S of the form (3.10). Substituting (3.10) into (4.12) and 
solving the resulting equation we obtain non-trivial forms of S for the following 
specific parametric restrictions (we omitted the uncoupled case 612 621 = 0), namely 

(i) ai=a 2 , 6 21 -6 11 (2-^)=0, (4.13) 

022 

(a) &21 = fon, 612 = &22, (4-14) 

/•••\ t 622611 ... 

(m) 621 = —r , (4-15) 

012 

61 o 

(iv) 01=03, 621-611(2-^)^0, (4.16) 

622 

and the respective S 1 forms are 

W (n)S = -f, 

(«<) 5 = -^, M^-4^±^). (4.17) 
6i 2 x' a;(6na; + 6121/) 

Now substituting the above forms of 5 into equation (3.3) and solving the resultant 
equation we obtain the corresponding forms of K. By making use of the ansatz 
mentioned in §3 6 we obtain following expressions for K, namely 

b 22 ( a l b 22 „ \+ 

(Hi) K = -x~^e ( ~^~~ a2)t , 

/• n „ (biix + bi 2 y)e b ^ x b ^ ,. 
(iv) ^=~r- \T> ( 4 - 18 ) 

(r-l) (fci2&2i+fcii(fci2-262 2 )) _ u _ & 
r (b 12 -b 22 )x ° llX ° 12 y 

with r is given in equation (4.5). It may be noted that the set (4.13)-(4.16) also 
includes the known time dependent integrable cases. 

Once R(= SK) and K are determined then one has to verify the compatibility 
of this solution with the extra constraint (2.8) which indeed gets satisfied in each 
one of the above four cases. Substituting the resultant integrating factors into (2.9) 
and evaluating the integrals we obtain the associated time dependent integrals of 
motion in the forms 

(■ \ t e~ a2t (bux + b 22 y) 

I= {a 2 + bll x + b22y y fl ^°' ^ 

db) f= 1 + (^ + W , fl2 ^ , (4.20) 
(bux + b 22 y) 

(ii) J= £ e («w-°i)*, (4.21) 

V 

(Hi) I = e ( ^~ a2) *a;~Sy, (4.22) 
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(iva) 

(ivb) 
(ivc) 



where 



di 

d 3 
e3 



(r - 1) d 3 x 
r d\ 



e 2 



^ ^ d\x 2 + e 2 d 1 



x {b 11 b 22 x + b\ 2 y - 6 12 ei) 



(r ~ 2)t ' 22 £2£L(r-2)t 

X b i2 e b i2 v ; 



/ = log 



2x *12 yd\ 



b 21 x a 2 di 
b\ 2 y b 12 



t, 



r ^ 0,2,(4.23) 
r = 2, (4.24) 



I = e "i2 x "i2 \b\ x d 2 x 2 + b\ 2 d x y 2 + b xx b X2 e±x , r = 0, (4.25) 



(612-622), ei = (621a; + b 22 y), d 2 = (6 i2 - 26 22 ), 

(612621 + b n d 2 ), e 2 = (bux + b 12 y), 

((bn-b 21 )x + 2d 1 y) 1 e 4 = {b 21 x + 2d lV ). (4.26) 



Finally we note that our method not only gives us a rather general set of inte- 
grate parametric relations (which includes all known cases) but also provides two 
independent integrals from which one can deduce the general solution for some 
cases. 



(iii) General solutions /Integr ability 

Interestingly, we find that for certain parametric choices we have two indepen- 
dent integrals (time independent as well as time dependent) and consequently one 
can express the general solution explicitly by using both of them. For example, 
let us consider the parametric choice given in (4.13) and the associated time de- 
pendent integrals given in (4.19) and (4.20). For this parametric choice, that is, 
a\ = a 2l b 2 \b 22 — b\\(2b 22 ~ b\ 2 ) — 0, we can also find the following time indepen- 
dent integral from (4.9), 



y 



{a x + b n x + b 22 y)' 



(4.27) 



Using the integrals / and I 2 , the general solution for the 2D LV system, (4.1), for 
the parametric choice (4.13) can be written as 



aiaie ait I 
x(t) = ^ 

buaigi + 622/2(51) b22 



and 



x(t) 



h(92) 622 - 611 92 



b 12 _ 



y(t) 



y(t) 



aie ait II 2 



(61101(51) b22 +b 22 I 2 gi) 



oi ^ (4.28) 



f£l2._ii 

U22 I 



622(61152 --Mfl^)" 22 ) 



-, Oi = 0, 



(4.29) 



where di = a^ 22 , 51 = (1 — Ie ait ) and g 2 = (t + I), respectively. Further, for 
the parametric choice given in (4.14), that is, 621 = 6n, 612 = 622, we obtain the 
general solution of the form 

u\ aia 2 e ait ai a 2 e a ^ 
x(t) = — r- , y(t) = — — , (4-30) 



{II-^)—.-gz(t) 



(//f a2 )^ -9s(t) 
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where gz{t) = {a\b 22 e a2t + a 2 b\\I\e ait ) . Depending on the signs and magnitudes of 
the system parameters oi, a 2 , &11, 612 and 622, the above solutions describe normal- 
ized interacting populations which asymptotically decay or grow or saturate. 

Similarly for all the other integrable cases identified in this section, one can 
derive the general solutions which are physically and mathematically relevant, often 
after some manipulations. The details will be presented elsewhere. 

(b) Application to non- autonomous system of first order ODE 

As we mentioned in the introduction, one of our motivations is to show that the 
procedure developed in §3 is also applicable to non-autonomous systems as well. 



(i) Example 1: Complex Riccati equation 

To demonstrate this point in brief let us consider the following first order non- 
autonomous equations 

x = ai(t)x - a 2 (t)y + /3i(t)(x 2 - y 2 ) - 2(3 2 (t)xy = </>i(t, x, y), 

y = ai (t)y + a 2 (t)x + [3 2 (t)(x 2 - y 2 ) + 2fa{t)xy = 4> 2 {t, x, y), (4.31) 

where di(t) and (ii{t) 1 i = 1,2, are arbitrary functions of i. Equation (4.31) describes 
the dynamics of two interacting species with time modulated parameters. This cou- 
pled equation is essentially the real form of the complex Riccati/Bernoulli equation 
z = a{t)z + j3{t)z 2 , where z = x + iy, a{t) = a 1 (t)+ia 2 (t) and 0(t) = 0i(i) + ilh(i). 
Substituting the form of 4>iS, i = 1,2, into the determining equation (3.2) and solv- 
ing the latter we obtain the following forms for S (generalizing the ansatz (3.10)), 
namely 

(y 2 - x 2 ) cos(w(i)) — 2xysin(w(£)) 

S\ = 



(x 2 — y 2 ) sa\{uo{t)) — 2xy cos(ui(t)) ' 
(y 2 - x 2 ) sin(w(i)) + 2xy cos(u(t)) 
(y 2 — x 2 ) cos(a>(t)) — 2xy sin(w(£)) ' 



(4.32) 



where Lu(t) = J a 2 (t)dt. Now inserting the above forms of Si and S 2 into equa- 
tion (3.3) and solving the resultant equations we obtain the following integrating 
factors, that is, 

e f ai (t)dt / \ 
K * = ( x 2 +y 2)2 [(y 2 - x2 ) sin(w(t)) + 2xy cos(a,(i))J , 

J ai (t)dt / \ 
^2 = {x 2 + y 2 ) 2 [i x - V ) sin(w(i)) + 2xy cos( W (t))J . (4.33) 

Substituting the complete sets (Ri(= SiKi), Ki), i = 1,2, into equation (2.9) and 
evaluating the integrals we obtain 

J ai (t)dt 

{x 2 + y 2 ) 
e f ai (t)dt 

7 2 = ( x 2 +y 2) ( xsm(uj(t))-ycos(uj(t)) ) +j 2 (t), (4.84) 



7l = 1 2 1 2\ \ xcos(w(t)) +ysm(w(t) ) +Ji(t), 
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where 

7iW= J (fa(t) cos(uj(t))-p 2 (t)sm(uj(t))yf a ^ dt dt, 

72(*) = y ^ 2 (t)cos(a ; (t))+/3 1 (t)sin(a ) (t))^ e J' Ql W dt *. (4.35) 

From the integrals I\ and 7 2 , the general solution for the equation (4.31) can be 
fixed easily in the form 

e f ai (t)dt / \ 
x(t) = ^(h - 7i W) cosMi)) + (h - 72 W) sin(w(t))J , 

f ai (t)dt / \ 
= 7s(t) (^(/i-7iW)sin(o;(t)) - (I 2 - 72(f)) cos(w(t))J , (4.36) 

where 73(f) = (/1 — 7i(f)) 2 + {h —72(f)) 2 - Again depending on the nature of the 
system parameters, the above solutions represent oscillatory or decaying or growing 
populations. 

(ii) Example 2: A predator-prey equation 

To demonstrate the theory for non-autonomous systems further, we consider 
another example which is a non-autonomous predator-prey equation, 

x = axx + j 1 e~ /3lt xy = fa, y = a 2 y + ^e'^xy = fa, (4.37) 

where on, Pi and 7,, i — 1,2, are arbitrary parameters. Substituting the form of 
fa's, i — 1,2, into the determining equation (3.2) and solving the latter we obtain 
the following form for S for the parametric choice 0i = (3 2 = 0, namely 

Si - J^-^ + ^y . (4.38) 
[ai — p)x + jie~ l:lt xy 

Now inserting the above form of Si into equation (3.3) and solving the resultant 
equation we obtain 

^K-fl+^e-^ 3 
V 

Substituting the forms R\{= S\K{) and K\ into equation (2.9) and evaluating the 
integrals we obtain 

h = (72a; - liv)^ t3t + {a 2 - fi) log x - {ax - (3) log y - (a 2 - atfjPt. (4.40) 

Unfortunately we could not find a second integral in this case. However, for the 
further parametric restriction, namely ct\ = a 2 = a we obtain S 2 of the form 

S 2 = (4.41) 

7i 

Substituting S 2 into equation (3.3) and solving the resultant equation we obtain 

K 2 = lie- at . (4.42) 
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Plugging the forms R 2 (= SK) and K 2 into equation (2.9) we obtain I 2 as 

/ 2 = (72^ - 7iJ/)e~ Qt . (4.43) 

From the integrals I\ and 7 2 , the general solution for the equation (4.37) with 
ai = a 2 = a and (3i = (3 2 = can be fixed easily in the form 

*(*) = bags ' y(t) = — «-S->>»> — • (4 - 44) 

(72-7ie 3=5 ) ( 72e s=73 - 7l ) 

It is obvious that depending upon the signs and magnitudes of the parameters a 
and (j the general solution cither decays or grows or saturates in the asymptotic 
limit. 

5. Prelle-Singer procedure for three coupled first order 

ODEs 

Next, we focus our attention on a system of three coupled first order ODEs of the 
form 

= P!(t,x,y,z) . = P 2 {t,x,y,z) = P 3 (t,x,y,z) 
Qi(t,x,y,z)' Q 2 (t,x,y,z)' Q 3 (t,x,y,z)' 

where Pi's and Qi's, i = 1,2, 3, are analytic functions in 2, y and 2; with coefficients 
in the field of complex numbers. Further, we assume that the ODE (5.1) admits 
a first integral I(t, x, y, z) — C, with C constant on the solutions so that the total 
differential becomes 

all = I t dt + I x dx + I y dy + I z dz = 0. (5.2) 
Now let us rewrite the equation (5.1) in the form 

P P P 

-±dt -dx = 0, -?-dt -dy = 0, -±dt - dz = 0. (5.3) 
Qi Q2 Qz 

Hence, on the solutions, the 1-forms (5.2) and (5.3) must be proportional. Multi- 
plying the first, second and third equations in (5.3) by the functions R(t,x,y,z), 
K(t,x,y,z) and M(t,x,y,z), respectively, which act as the integrating factors of 
the corresponding equations, we have on the solutions that 

dl = (J20i + K(j) 2 + M4> 3 )dt - Rdx - Kdy - Mdz = 0, (5.4) 

where = Pi/Qi, i = 1,2, 3. Comparing equations (5.4) and (5.2) we have, on the 
solutions, the relations 

It = {Rfa + K<j> 2 + M0 3 ), I x - -R, I y - -K, I z - —M. (5.5) 

The compatibility conditions between the equations (5.5) provide us the follow- 
ing determining equations for the integrating factors R, K and M: 

Rt + (f>iR x + <j) 2 R y + foR z = -{R<pi x + K<p 2x + M<j> 3x ), (5.6) 

K t + <t>iK x + faKy + (f> 3 K z = -(R<f>iy + K4> 2y + M4> 3y ), (5.7) 

M t + 4>iM x + <j> 2 M y + 4> 3 M Z = -(R4> lx + Kcf> 2z + M(f> 3z ), (5.8) 

R y = K x , R Z = M X1 K z = My. (5.9) 
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On the other hand integrating equations (5.5), we obtain the integral of motion, 

d 



I = n + r 2 + r 3 



-I 



where 



r\ = 



r 3 = 



R<p!+ K(f> 2 + M<p 3 )dt, r 2 



M + — [ n + r 2 + r 3 



^jdt, r * = -j 



dz, 



(5.10) 



K+—{r 1 +r 2 ) )dy. 



Naturally, for the complete intcgrability of Eq. (5.1) we require three independent 
integrals and so three independent sets of integrating factors (Ri,Ki,Mi), i — 



1,2,3. 



6. Method of Solving Determining Equations 

The determining equations (5.6)-(5.9) are more complicated than the two dimen- 
sional case discussed in §3 and so to simplify the determining equations we introduce 
the transformations 

R = SM and K = UM, (6.1) 
where S and U are functions of t, x, y and z, so that the equations (5.6)-(5.9) become 



D[S] 


= S{S<p lz + Ufo z + 4 


>3z) - {S<j>l x + U(f> 2x + (f>3 x ), 


(6.2) 


D[U] 


= U(S<p lz + Uc/> 2z + ^ 


^z) - {S<f>ly + U(f> 2y + <j) 3y ), 


(6.3) 


D[M\ 


= -M(S<f> lz + U<j> 2z - 


f <hz), 


(6.4) 


M x 


= SM Z + MS Z , My 


= UM Z + MU Z , 


(6.5) 


U X — Sy 


= SU Z -US Z , D = 


d d d d 

-tt. + <Pi-K- + <P2Tr + 03^-- 
at ox oy Oz 


(6.6) 



One may note that two of the determining equations are still in coupled form and 
the transformations are natural extensions of the two coupled case. Of course, one 
may also consider alternate possibilities, that is, either R = SK and M — UK or 
K = SR and M = UR. However, such possibilities again lead to same results. 

In the two coupled case (2.1), we divided our analysis into two categories, vide 
§3 a (time- independent integrals) and §3 6 (time-dependent integrals). However, in 
the present case we divide our analysis into five categories, that is, (i) I x = and 
It, ly, h + 0, (ii) I y = and I t , I x , I z ^ 0, (iii) I z = and I t , J x , I y ^ 0, (iv) 
It = and I Xl I y , I z and (v) I t , I x , I y , I z 7^ 0. We intend to proceed in 
this way because we observed the absence of a dynamical variable in some integrals 
in certain specific dynamical systems of the type (5.1). We try to identify these 
cases first. In fact, proceeding in this way we are able to formulate a condition 
on the system variables and if the given system satisfies this condition one can 
conclude that the given dynamical system has the integral without that respective 
variable. Since we have four variables, t,x,y and z we consider each one of the cases 
separately and treat none of the variables being absent as the fifth independent case. 
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Since we are dealing with a system of coupled three first order ODEs, the com- 
plete integrability is guaranteed by the presence of 2 time-independent integrals 
(whereupon the system can be reduced to a single quadrature) or 3 time dependent 
ones (in which case the solution can be obtained in an algebraic way (Bountis et 
al. 1984)). In the following we will search for such integrals. 

(a) Caes 1: I x = and I t , I y , I z ^ 

In the case I x = 0, we have R — (vide equation (5.5)) which in turn implies 
that (i) either 5 = and M ^ or (ii) 5 ^ and M = as can be seen from (6.1). 
In the former case, 5 = 0, M ^ 0, one can easily fix the form of U, from equations 
(6.2) and (6.6) as 

U = -^, U x = 0. (6.7) 

02a; 

On the other hand the choice M = and 5^0 when R = leads to the case 
where one of the dynamical variables become uncoupled (see equations (5.6)-(5.9)) 
and so effectively a system of two coupled first order ODEs results in, which we 
have already discussed. So this choice is not considered further. Inserting the above 
form (6.7) into (6.3) we arrive at the condition 

03a-(02a;i + fofoxy + 0302xz ~ 02z03a- ~ 02x022/) ~ 02a;(03a;i 

+0203xi/ + 0303xz - 02x03j/ - 03x03z) = 0. (6.8) 

The condition (6.8) gives us the integrable cases for which the system possesses the 
integrals of motion with I x = 0. Now substituting (6.7) into (6.4) we obtain the 
following determining equation for M, that is, 

D[M] = M(^0 2z -0 3z ). (6.9) 

02a; 

Again to solve equation (6.9) one has to make a suitable ansatz for M. Choosing 
appropriate ansatz for M and solving the equation (6.9) one can get explicit form 
for M. Once M is known, the integrating factors can be fixed from the relations 
K = UM and R = SM = 0. Now plugging the forms of R, K and M into Eq. (5.10) 
and evaluating the integrals one can construct the integrals of motion for the given 
system. 

(b) Case 2: I y = and I t , I x , I z ^ 

The determining equations and conditions can be fixed in a similar manner for 
this case, I y = and I t , I x , I z ^ 0, with the replacement (5, U, M, 0i, 02, 03, x, y, z) 
— ► (U, 5, M, 02, 0i, 03, y, x, z) in the above analysis. 

(c) Case 3: I z = and I t , I x , I y ^ 

In the present case with the form of the integrating factors R = SM and 
K = UM, I z = implies M = and so R — and K = as well, leading 
to an integral of motion which turns out to be constant. Therefore in this case 
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we consider the other possibility R — SK and M = UK and proceed as above. 
The final results are obtained with the replacement (S, U, M, 0i, 02, 03, x, y, z) — ► 
(U,S,K, <t>3,<f>i,<h^,x,y) in case 1. 

(d) Case 4: It = and I x , I y , I z ^ 
Next, in the time independent case I t = the first equation in (5.5) gives 

s * -te + W. (6 . 10) 

Substituting this form of 5 into (6.3) and (6.4) we get the following form of deter- 
mining equations for U and M, 

D[U] = ^ ^^ {Ufay - <j> lz ) + U{U<hz + foz ~ <j>2y) ~ foy, (6.11) 

D[M] - M ^ + ^ fa-faU-faz). (6.12) 

To solve the equations (6.11) and (6.12) we adopt the following methodology. 
To start with, in order to solve (6.11), we consider U in the form 

u= Mx,y) + B 1 (x,y)z (613) 

A 2 (x,y) + B 2 (x,y)z 1 

where A^s and Bj's, i = 1,2, are arbitrary functions of x and y. Substituting (6.13) 
into (6.11) and equating the coefficients of different powers of z to zero, we get a 
set of determining equations for the functions A^s and Si's, i = 1, 2,. Solving these 
determining equations we obtain explicit expressions of the functions A^s and B^s, 
i = 1,2, and consequently the associated function U . 

Now substituting the forms of U into equation (6.12) and solving the resultant 
equation we obtain the corresponding form of M. To solve the determining equation 
for M we again seek the ansatz of the form M = -r-r-, — ^-^7 — r-vF , where Ud is the 
denominator of U. Once U and M are fixed then one has to verify the compatibility 
of this set (S, U, M) with the constraint equations (6.5)-(6.6). Now substituting R(= 
SK), K(= UM) and M's into equation (5.10) one can construct the associated 
integrals. Finally, one can proceed with a more generalized ansatz than (6.13) if 
need arises. 



(e) Case 5: It, I x , Iy, Iz case 

Solving the determining equations (6.2)-(6.4) are naturally more tedious with 
none of the variables (t, x,y, z) absent in I, when compared to the earlier cases. 
To start with one may use the following simple ansatz to solve the determining 
equations (5.6)-(5.8), that is, 

R = Ai(t,x,y) +B 1 (t,x, y)z, K = A 2 (t, x, y) + B 2 (t, x, y)z, 
M = A 3 {t,x,y) + B 3 (t,x,y)z, (6.14) 

where A^s and Bj's i — 1,2,3, are arbitrary functions of t, x and y. Depending 
on the nature of the equation (5.1), one may work with more general forms like a 
rational one. 
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7. Three coupled ODEs - Applications 

(a) Example: Rossler system 
Let us consider the Rossler system (Rossler 1976) 
dx dy dz 

— = -y-z = (f> 1 , — = x + axy = <f) 2 , — = a 2 + xz + a 3 z = (j> 3 , (7.1) 
at at at 

where a^'s, i = 1,2,3, are arbitrary parameters. Several works have been devoted 
to study the dynamics of this equation. Very recently, Llibre & Zhang (2002) and 
Zhang (2004) have studied equation (7.1) using the so called Darboux method and 
obtained conditions for integrability. We now apply our above method to system 
(7.1) and explore new integrals, if they exist. 

(i) I x = and I t , I y , I z ^ 
Substituting (7.1) into (6.8) we get 

zai + a 2 = (7.2) 

From equation (7.2) we conclude that ct\ = a 2 = so that from (6.7) we get 
U = —z. The determining equation for M turns out to be 

M t + xM y + z(x + a 3 )M z = -M(x + a 3 ), (7.3) 

in which we have taken M x = (since I x = 0). A simple solution for (7.3) is 
M = 1/z. Making use of the explicit forms of U and M and with the parametric 
restriction a\ — a 2 = we conclude that i? = 0, K = — 1, M = -. Now substitut- 
ing the functions R, K and M into equation (5.10) and evaluating the integrals we 
obtain the following integral of motion 

I = y + a 3 t-\og(z). (7.4) 

The integral (7.4) with a 3 = has already been given by Llibre & Zhang (2002) 
and Zhang (2004). The integral (7.4) with a 3 ^ is new to the literature at least 
to our knowledge. 

(ii) I y = 0, I t , 4, I z ^0 and I z = 0, I t , I x , I y ^ 

Proceeding appropriately we could not find any integrable case in the Rossler 
system belonging to these categories. 

(hi) I t = and I x , I y , I z ^ 

In this case the function S can be fixed in the form (vide equation (6.10)) 

_ a 2 + (x + a 3 )z + (x + aiy)U 
y + z 

Substituting (7.1) into (6.11) we get 

nrrn a 2 + (x + a 3 )z + (x + ary)U , ur 

D\U\ = (U - 1) + (x + a 3 - ai)U. (7.6) 

y + z 
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Substituting (6.13) into (7.6) and solving the resultant equations wc obtain non 
trivial forms of U for the specific parametric restriction ot\ = a 2 = a 3 = and then 
making use of the U forms into (7.5) we obtain (Si, U\) = (x, y), (S2, U 2 ) = (0, —z). 
Now substituting the forms of U into equation (6.12) and solving the resultant 
equation we obtain Mi = —1 and M 2 = —e~ v . Now inserting the integrating 
factors i?i's(= SjM,), Ki's(= UiMi) and Mi's, i = 1,2, into (5.10) and evaluating 
the integrals we arrive at the expressions 

h=x 2 +y 2 + 2z, I 2 =ze-y. (7.7) 

These two integrals are already known (Llibre & Zhang 2002; Zhang 2004). 



(W) I t , I X , Iy, I Z ^0 

In this case, substituting <^'s, i — 1,2,3, into (5.6)-(5.8) and solving the resul- 
tant system of equations with the ansatz (6.14) we obtain the integrating factors 
for the parametric choice cti = a 3 = 0, that is, R = —x, K = —y, M = —1 and 
the corresponding integral of motion takes the form 

I = (x 2 + y 2 + 2z -2a 2 t). (7.8) 

The integral (7.8) has also been reported by Zhang (2004). 

We conclude this section by mentioning that our studies reveal that the system 
(7.1) possesses a time dependent integral for the parametric choice a 3 ^ 0. 

(6) Example:2 3D-LV system 

Let us consider a 3D LV model for competition between three populations whose 
dynamical evolution is determined by the following equations (Cairo 2000) 

x = x(ai + aix + hy + c\z), y = y(a 2 + a 2 x + b 2 y + c 2 z), 

z = z(a 3 + a 3 x + b 3 y + c 3 z) (7.9) 

where a.i, a^, bi and Cj, z = 1,2, 3, are arbitrary parameters. Needless to say 3D-LV 
system is one of the challenging problems and a testing ground for several analytical 
methods. In the following we identify the integrals of motion for certain specific 
parametric choices in (7.9) using the procedure given above. For the purpose of 
demonstration, in the following we present our results only for a couple of cases 
and a detailed analysis will be presented separately. 



(i) case 1: 

For the parametric choice at — a, bi — b and a — c, i = 1,2,3, we find the 
following three complete sets of integrating factors (B4, Ki,Mi), i = 1, 2, 3, 

e (a 3 -a 2 )t z p (a 3 -a 2 )t 

Ri=0, K 1 = - , Mi = , 

V V 

e (a 3 -a 1 )t z e (a 3 -cti)t 



R 2 = 5 , K 2 = 0, M 2 = 

x z x 



e ait h(v z) p" 1 * 
R 3 = — T' K 3 = a 3 — , M 3 = a 2 — , (7.10) 

x z ax ax 



Article submitted to Royal Society 



Integrability and Linearization 



19 



where h(y,z) = (a 2 a 3 + ba 3 y + ca 2 z). Now substituting the functions i?, K and 
M into equation (5.10) one can obtain the following integrals of motion 

g(<*2— Oi3)t z e (a 1 -a 3 )t z 

h = , h = , 

V x 

e ait {aia 2 a 3 + a 2 a 3 ax + ct x a 3 by + a\a 2 cz) 
i 3 — . I'- 11 ) 

x 

From the integrals h,I 2 and I3, we can deduce the general solution for the equation 
(7.9) for the parametric choice a» = a, bi — b and a = c, i = 1, 2, 3, of the form 



z(t) - 

y(*) = 

z(t) = 



aia2Ct3-fie Ql * 




— {aa 2 a 3 he ait + ba\a 3 I 2 e a2t + ca\a 2 hl 2 e a3t ) 
a\a 2 a 3 I 2 e a2t 


/3/1 


- (aa 2 a 3 I 1 e ait + botia 3 I 2 e a2t + coi\a 2 I\l 2 e ait ) 
etia 2 a 3 lil 2 e a3t 


/3/1 


— (aa 2 a 3 I\e ait + ba\a 3 I 2 e a2t + cotia 2 I\I 2 e a3t ) 



(7.12) 



For the above choice of parameter the system (7.9) is a completely integrable one. 



(ii) case 2: 

For the parametric choice on — a, i — 1,2,3, 03/3 = —a 2 /2 = 01, ci = 
— C3, bi = —63 and c 2 = b 2 = we obtain two complete sets of integrating factors 
of the form 

Ri = e~ iat V 2 z 1 K x = 2e- iat xzy, M 1 = e- iat xy 2 , 

R 2 = ai ay{%a Y x - 3c lZ + 2b 1 (y - — z))e" 3at , 

a 

4 

K 2 = a(3ajx 2 + bjy 2 + a lX (4b iy - Cl z(3 + -b iy ))e- 3a \ 

M 2 = -a lCl xy(3a + 2b iy ) e - 3at . (7.13) 

Now substituting the integrating factors R 7 K and M into equation (5.10) one can 
obtain the corresponding integrals of motions, namely 

h = e- 4at xy 2 z, 

I 2 = e- 3at ((auHxyiSaz - 3a lX - 2b lV + ^hayz) - \b\y^j • (7.14) 

We could not find the third integrating factor (R 3 , K 3 , M 3 ) within our ansatz, 
(6.14), and further detailed exploration is needed to conclude whether it exists or 
not to find the third integral if it exists. 



8. Prelle-Singer procedure for n (>3) coupled first order 

ODEs 

The above procedure to find integrating factors and integrals can be extended in 
principle to a system of n coupled first order ODEs (n > 3). In this cases, we 
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get n determining equations for the n integrating factors along with n(n — l)/2 
constraint equations, which one can solve algorithmically by following the above 
procedure. Further, the same procedure can also be applied in principle to any 
higher order as well as coupled higher order equations. This is because, any higher 
order equation can always be rewritten equivalently as a system of first order ODEs. 
For example, the second order equation x = f(x,x,t) can be written in the first 
order form as x — y, y — f(y,x,t). Then the determining equations for this first 
order form, namely equations (3.2)-(3.4) can be related to the determining equations 
(2.6)-(2.8) given by us earlier (Chandrasekar et al. 2005) for the second order form 
x = f(x,x,t). However, one has to note that each one of the procedures has its 
own merits and demerits, for example in constructing nonstandard Hamiltonian 
structures. Similarly, a straightforward first order form for the third order equation 
x = f(x, x, x, t) is x = y, y = z, z = f(x, y, z, t). In this case also the determining 
equations for the first order form can be related to that of the third order ODEs as 
given by us earlier (Chandrasekar et al. 2006). Again analyzing third order equations 
as such has its own practical advantages. This analogy can also be extended to 
coupled higher order equations as well in principle, though in the actual analysis one 
may have to make a judicious choice of which one of the methods is advantageous 
for investigating the integrability aspects. 



9. Linearization 

In this section, we describe a procedure to deduce the linearizing transformations 
from the known integrals and illustrate the theory with an example. 
Let us assume that equation (2.1) admits the following integral 

I = F(x,y,t). (9.1) 

Now let us split the function Fi in the form, 

Now we identify the function G\ as the new dependent variable and the integral of 
G2 over time as the new independent variable, that is, 



w = G 1 (t,x,y), r = [ G 2 (t',x,y)dt', (9.3) 

J o 

We note here that the integration on the right hand side of (9.3) leading to r 
can be performed provided the function Gi is an exact derivative of t, that is, 
G2 = 4;T(t, x, y) = xt x + yTy + r t . In terms of the new variables, equation (9.2) can 
be modified to the form 

'-«(£)• ,M) 

Inverting the relation (9.4) suitably, one can obtain a linear equation, 

dw t 

- = I, (9.5) 
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where I is a constant. Or equivalcntly 

dw t du . . 

-=u = I, -=0. (9.6) 

Equation (9.6) is the corresponding linear equation of (2.1). From equations (9.3) 
and (9.6) we have the following linearizing transformation for (2.1), namely 

w = G 1 (t,x,y), u = F(x,y,t), r= f G 2 (t' ,x,y)dt' . (9.7) 

J o 

In this case the new variables u>, u and r helps us to transform the given system of 
coupled first order nonlinear ODEs into a linear system of coupled first order ODEs 
which in turn leads to the solution by trivial integration. The above procedure can 
also be extended to more than two coupled first order ODEs straightforwardly and 
we do not present the details here. 



(a) Example 1: 

To illustrate the underlying ideas let us consider the two dimensional LV system 
(4.1) with the specific parametric choice, &n = 621 and &12 = 622, namely 

x = x(ai + bux + b 22 y), y = y{a-2 + bux + b 2 2y)- (9.8) 
Let us consider the following first integral for equation (9.8), namely 

j = V e (a 1 -a 2 )t_ /gg) 
X 

Now rewriting equation (9.9) using (9.8) in the form (9.2) we get 



022 dt 



[(i- 






ai J 



(9.10) 



Then 



w= (I + ^lW, T =- b Jl e ^. (9.11) 
\x ai / a 2 

From equations (9.11), (9.9) and (9.6) we have the following linearizing transfor- 
mation for (9.8), namely 

hl\ P *it „ - l P (ai-a 2 )t ^__^22 a 2 t 



w= - + — e ait , u= ± e ^- a2 >\ T =-—e a >\ (9.12) 
\x ai / x a 2 

In this case the new variables w, u and r helps us to transform the given system of 
coupled first order nonlinear ODE, (9.8), into a linear system of coupled first order 
ODEs of the form (9.6). The general solution can then be deduced straightforwardly. 
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(b) Example 2: 

Similarly, for the specific parametric choice a\ = a 2 = 0, 612 = — &22 and 
&21 = 3&n, equation (4.1), that is, 

x = x(bux + b 12 y), y = y(3bux - h 2 y), (9.13) 

can be transformed to linear equation of the form (9.6) by the following linearizing 
transformations, namely (for the integral see equation (4.20)) 

1,2 (bm + bux) 1 

W =o t -^7 77 7 72' U = t+ I 1 ' T = t - ( 9 ' 14 ) 

2 2bux(bi2y - biix) 2 b\\x-b\iy 



10. Conclusion 

In this paper, we have modified the Prcllc-Singer procedure such that it is ap- 
plicable to both autonomous as well as non-autonomous system of coupled first 
order ODEs. We have also developed systematic procedures for finding both time 
independent and time dependent integrals for them. From this analysis we have 
answered the following open questions (i) How can the PS method be extended 
to non-autonomous system of coupled first order ODEs? (ii) How to find the sec- 
ond or the time dependent integrals for the given coupled first order ODEs and 
(iii) How can this procedure be generalized to higher dimensions in order to find 
first integrals? We have also shown that the determining equations for the time 
independent integral in the two coupled equations of our method coincides with 
the determining equations derived by Prelle and Singer in their original paper. We 
have illustrated this procedure with physically interesting examples, namely two 
dimensional LV system, Rossler system and 3D-LV system and identified several 
integrable cases. Further, we have developed a linearization procedure for coupled 
first order ODEs. Finally, we note that the procedures which we have developed in 
this paper, namely both the extended PS procedure and linearization can also be 
extended to any number of coupled first order ODEs. 

The work of MS forms part of a research project sponsored by National Board 
for Higher Mathematics, Government of India. The work of ML forms part of a 
Department of Science and Technology, Government of India sponsored research 
project and is supported by a DST Ramanna Fellowship. 
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